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Abstract 

The  basic  equations  and  some  fundamental  concepts  of  mag¬ 
netoaerodynamics  are  discussed.  As  an  application,  the  problem  of  the 
flow  near  the  stagnation  point  of  a  body  of  revolution  is  reviewed  and  an 
exact  solution  is  given.  Data  are  presented  for  the  heat  transfer  coef¬ 
ficient  at  the  stagnation  point,  and  for  its  gradient  in  the  stagnation  point 
region.  The  magnetic  field  strength  required  to  accomplish  an  appre¬ 
ciable  reduction  of  aerodynamic  heating  in  hypersonic  flight  is  discussed, 
for  the  case  in  which  ionization  is  due  to  thermal  motion.  Alternatively, 
methods  which  involve  electrical  breakdown  of  the  air  are  considered. 

Introduction 

At  hypersonic  speeds,  exceeding  a  flight  Mach  number  of 
approximately  15,  the  air  behind  a  normal  shock  reaches  a  temperature, 
which  is  sufficient  to  ionize  the  air  to  an  appreciable  extent.  It  is  then 
possible  to  alter  the  characteristics  of  the  flow,  by  means  of  externally 
applied  magnetic  fields,  (1)  to  (8).  Due  to  the  motion  of  the  ionized 
gas  (plasma)  through  the  magnetic  lines  of  force,  currents  are  induced 


*)  Member  of  Senior  Staff 

**)  Number  in  parentheses  indicate  References  at  end  of  paper. 
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in  the  gas,  which  in  turn  give  rise  to  a  force  j  x  B,  per  unit  volume, 

— >  — ► 

acting  on  the  gas.  In  this  expression,  j  is  the  current  density  and  B 

the  magnetic  induction.  The  plasma  currents  in  turn  give  riee  to  a 

— ^ 

magnetic  field,  which  adds  to  the  field  produced  by  external  currents.  B 
represents  the  total  field.  The  theoretical  analysis  of  the  ensuing  flow 
therefore  requires  in  general  a  solution  of  Maxwell's  equations  for  the 
electromagnetic  field,  coupled  with  the  equations  of  fluid -mechanics  con¬ 
taining  the  additional  force  term.  The  equation  of  energy  is  also  modi¬ 
fied,  in  a  manner  which  will  be  discussed  below. 

Confining  the  present  considerations  to  the  stationary  flow 
of  a  fluid  of  constant  properties,  the  usual  equations  of  magnetohydro¬ 
dynamics  take  the  following  form: 

curl  H  =  j 


curl  E  =  0 


(1) 
(2) 

div  B>  =  0  (3) 

_ v  _ 

where  H  and  E  are  the  magnetic  and  eJectric  field  strengths,  respec-r 

tively.  Maxwells  equations  are  supplemented  by  the  constitutive  equation 

— *  — 

B  =  pH  (4) 

and  Ohm's  law  for  moving,  isotropic  media, 

— >  — »  — »  — * 

j  =  <r  (E  +  v  x  B)  (5) 

where  p  is  the  magnetic  permeability,  <r  the  conductivity  and  v  the 
velocity.  Furthermore, 

div  v  =  0  (6) 

since  a  fluid  of  constant  density  is  assumed.  Supplementing  the  Navier- 
Stokes  equation  by  the  magnetic  force  term, 


Dv* 


pI5f 


,  2  —>  ->  — 
=  -  grad  p+T/V  v  +  jxB 


(7) 


(p  =  density,  p  =  pressure,  =  viscosity) 
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Eqs.  (1)  to  (7)  suffice  to  determine  for  instance  the  velocity 
and  the  current  distribution  in  the  fluid.  For  the  purpose  of  computing  the 
distribution  of  temperature  and  rates  of  heat  transfer,  consideration  must 
be  given  to  the  equation  of  energy  ' , 


pc§^  -  .v  \>  2  T 


1  ~>2 

7  > 


♦  I 


(8) 


(c  =  specific  heat  per  unit  mass ,  T  =  temperature,  =  thermal 
conductivity),  where 


?v. 

i 

Tij 


T.  . 

y 


Jx. 


(9) 


•  ^  ♦  *  < 
The  terms  on  the  left  of  Eq.  (8)  represent  the  rate  of  heat  removed  per 

unit  volume  by  convection  and  by  heat  conduction,  respectively.  The  terms 

on  the  right  are  the  rates  of  heat  added  through  Joule  heating  and  through 

viscous  dissipation. 

Stagnation  Point  Flow 


As  an  illustration,  the  rotationally  symmetric  stagnation 
point  flow  is  considered.  The  corresponding  two-dimensional  case  was 
treated  by  Neuringer  and  Mcllroy  (10).  In  the  present  analysis,  a  ro¬ 
tationally  symmetric  magnetic  field  is  assumed,  which  vanishes  at  a 
large  distance  forward  of  the  stagnation  point  and  has  a  normal  component 
on  the  surface  which  is  constant.  It  follows  from  simple  symmetry 
arguments  that  a  solution  exists,  for  which 


*)  A  full  discussion  of  the  equation  of  energy  in  magnetohydrodynamics, 
applicable  to  more  general  conditions  than  needed  here,  is  due  to 
Goldstein  and  Chu  (9). 


Page  4 


j 


where  r ,  z  and  </>  indicate  the  corresponding  components  in  cylindrical 


coordinates.  It  is  seen  that  the  current  flows  in  circular  loops  about  the 

— ^  — * 

z-axis.  The  magnetic  force  j  x  B  is  in  the  meridional  plane,  and  its 
direction  is  such  as  to  impede  the  motion  of  the  fluid. 

The  trial  solution 
vr  =  r  V(z) 
vz  =  -2  i(z) 

- - B.  =  r  » '(*)  -  -  -  -  - - 


Bz  =  -  2  g(z) 

—  2 
p  -  p  =  h(z)  +  r  k(z) 

T  -  T  =  l{z)  +  r^m(z) 
where  f,g,h,k,l  and  m  are  functions  of  the  coordinate  z  alone,  and 
where  p  and  T  are  the  pressure  and  temperature  at  the  stagnation 
point,  satisfies  Eqs.  (1)  to  (9)  identically,  resulting  in  a  set  of  differen¬ 
tial  equations  for  these  functions.  Full  details  are  given  elsewhere  (11). 

The  following  boundary  conditions  are  imposed  on  the  sur¬ 
face,  i.e.  at  z  =  0, 


(11) 


(a)  the  normal  component  of  the  total  magnetic  field  is 
constant,  i.e.  B  (0)  =  B  where  B  is  a  given  con- 
stant. 

(b)  the  velocity  vanishes,  i.e.  v  (0)  =  v  (0)  =  0. 

r  z 

(c)  the  temperature  of  the  surface  is  constant,  i.  e. 

T  (0)  =  T. 
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At  a  large  distance  from  the  surface,  the  following  con¬ 
ditions  are  imposed, 

(d)  the  magnetic  field  vanishes.  In  particular,  it  is 
required  that  B^oo)  =  0. 

(e)  the  velocity  and  pressure  distributions  approach  the 

corresponding  distributions  for  an  inviscous,  non¬ 
magnetic  stagnation-point  flow.  In  particular, 
vr(<>0)  =  ar*  where  the  constant  "a"  is  determined 

by  the  characteristics  of  the  flow  outside  the  stagna¬ 
tion  point  region. 

(f)  the  temperature  is  independent  of  the  radius.  In 
particular,  T(  <)  =  T  where  T  is  a  given  constant. 

Numerical  results  were  obtained  by  means  of  a  digital  com¬ 
puter.  Introducing  the  dimensionless  quantitier 


£  =  \  p<ra  s 

(12) 

f(£)  =v 

(13) 

G<£)  =  -  |  g(x) 

(H) 

typical  distributions  of  velocity  and  magnetic  induction  are  plotted 

These  results  depend  upon  two  dimensionless  parameters,  vi*. 

in  Fig.  1. 

*  =  a(B'2/2ji)’1 

(15) 

P  =  por  ^  /  p 

* 

(16) 

In  aerodynamic  applications,  ^  and  p  are  invariably  very 
small,  typically  of  the  order  of  10"^.  The  ratio 
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y  =  pn*  =  ^  <17> 

which  is  a  measure  for  the  ratio  of  the  magnetic  forces  to  the  inertial 
forces,  however  is  generally  of  order  one.  Letting  •(  and  {3  tend  to 
zero  gives  rise  to  a  boundary  layer  problem  (11).  The  velocity  and 
the  magnetic  field  distributions  then  depend  only  upon  the  single  para¬ 
meter  y  .  Figs.  2  and  3  represent  the  two  velocity  components  out¬ 
side  this  boundary  layer  =  F).  Similarly,  Figs.  4 

and  5  give  the  velocity  distribution  in  the  boundary  layer,  where 
C2  =  P'1/2£and  F2  =  P"1/2F. 

As  one  recognizes  from  the  velocity  distributions,  an 
important  effect  of  the  magnetic  field  consists  in  displacing  the  .low 
in  the  direction  of  the  positive  z-axis.  The  thickness  of  this  magnetic 
cushion,  defined  by 

l  -  Z  ;  Z  =  lim  (z  -  f(z)/a)  (18) 

z-xx. 

is  plotted  in  Fig.  7. 

Defining  the  Nusselt  number 

Nu(r)  =  h(r)  i  (i)1/Z  (19) 

where  h(r)  is  the  heat  transfer  coefficient,  one  obtains 

Nu(r)  =  Nu(0>  +  Nu^  — r - 

c(T-T) 


(20) 
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The  coefficient*  Nu^  and  Nu^  have  been  computed  for  a  Prandtl  number 
of  0.70  and  are  given  in  Fig.  8.  With  increasing  magnetic  field  strength, 
the  rate  of  heat  transfer  at  the  stagnation  point  is  seen  to  decrease,  a  re¬ 
sult,  which  is  in  qualitative  agreement  with  the  computations  by  Neuringer 
and  Mcllroy  (10)  for  the  two-dimensional  case.  In  the  non --magnetic  case, 
which  corresponds  to  v  =0,  the  result  agrees  with  data  computed  by 
Cohen  and  Reshotko  (12). 

Nu^  increases  with  increasing  field  strength,  indicating 
an  increased  slope  of  the  heat  transfer  rate  considered  as  a  function  of 
the  radius.  In  applications  to  aerodynamic  heating,  this  is  a  detrimental 
effect.  The  cross-over  point,  at  which  the  advantage  of  a  reduced  heat 
transfer  at  the  stagnation  point  is  cancelled  by  an  increased  heat  transfer 
at  some  distance  from  this  point,  depend*  upon  the  details  of  the  applica¬ 
tion. 

Required  Magnetic  Field  Strength 

It  is  seen  that  -y  must  be  large  (of  the  order  of  100),  if 
an  appreciable  reduction  in  heat  transfer  at  the  stagnation  point  is  to  be 
achieved.  The  same  does  not  need  to  be  true,  however,  for  configura¬ 
tions  of  the  magnetic  field  differing  from  the  one  presently  considered. 

The  contribution  of  the  plasma  currents  to  the  total  magnetic  field  is  sub¬ 
stantial,  as  seen  from  the  curves  in  Fig.  6  where  the  magnetic  field  out¬ 
side  the  boundary  layer  is  plotted  in  terms  of  G^(£j)  s  G(£  ). 

The  required  field  strength  depends  critically  upon  the 
conductivity  of  the  medium.  The  conductivity  of  air  in  thermodynamic 
equilibrium  has  been  estimated  for  temperatures  up  to  24,  000°K  (Fig.  9). 

The  quotient  p/pQ  in  this  figure  is  the  ratio  of  the  density  to  the  density 
at  n.t  p.  These  data  have  been  computed  from  the  electron  collision 
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cross  sections  (13)  and  the  calculated  equilibrium  composition  of  air  (14). 
Details  are  given  in  an  Appendix  to  this  paper.  It  must  be  emphasized 
that  in  view  of  the  uncertainty  in  the  presently  available  data  on  the  appro¬ 
priate  cross  sections,  the  degree  of  accuracy  is  small. 

In  order  to  permit  a  quick  orientation,  the  conductivity  of  air 
behind  a  normal  shock  is  plotted  in  Fig.  10  directly  as  a  function  of  the 
flight  velocity  V  and  of  the  altitude  H  above  sea  level.  In  Fig.  11,  the 
magnitude  of  the  magnetic  induction  B  is  plotted,  which  is  required  to 
make  the  ratio  of  the  magnetic  forces  to  the  inertial  forces,  i.e.  aB^l/pu 
equal  to  unity.  Here,  a,  p  and  u  are  the  conductivity,  density,  and 
velocity  behind  a  normal  shock  in  air  (taking  into  account  dissociation  and 
ionization).  1  is  a  characteristic  length  of  the  aerodynamic  body,  and  was 
taken  as  lm.  With  these  assumptions,  a  definite  field  strength  is  associa¬ 
ted  with  any  given  velocity  and  altitude.  The  curve  labeled  by  the  equation 
u  =  V  is  discussed  in  the  Appendix. 

It  appears  from  these  data,  that  large  field  strengths  are 
required,  if  the  heat  transfer  at  the  stagnation  point  is  to  be  reduced  sub¬ 
stantially,  at  altitudes  for  which  aerodynamic  heating  is  critical.  Differ¬ 
ent  configurations  of  the  magnetic  field  may  be  more  advantageous. 

In  some  cases,  the  temperature  of  the  electron  gas  may  be 
many  times  larger  than  the  temperature  of  the  gas,  corresponding  to  an 
effective  increase  of  the  conductivity.  This  is  the  case  for  instance  with 
microwave  breakdown  and  glow  discharges.  The  device  shown  in  Fig.  12 
contains  a  large  number  of  electrodes  of  alternating  polarity,  which  are 
imbedded  in  an  aerodynamic  surface.  It  has  been  found  possible  to  main¬ 
tain  a  glow  discharge  in  a  supersonic  airstream,  ionizing  thereby  the  air 
in  the  boundary  layer.  A  typical  current  voltage  characteristic  is  given  in 
Fig.  13. 
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Appendix 

The  electrical  conductivity  of  air  in  thermodynamic  equilib¬ 
rium  at  temperatures  up  tc  about  7,000°K  is  known  with  considerable  re¬ 
liability  from  shock  tube  experiments  reported  by  Lamb  and  Lin  (13).  At 
these  temperatures,  the  contribution  of  the  positive  ions  to  the  total  electron 
collision  probability  is  still  relatively  small,  but  becomes  increasingly 
important  at  higher  temperatures. 

It  appears  that  no  detailed  experimental  data  are  available 
at  present  for  temperatures  in  excess  of  7,000°K.  On  the  other  hand, 
many  of  the  contemplated  applications  of  magnetoaerodynamics  involve 
consider, '•<'  higher  temperatures.  In  this  paper,  estimates  of  the 
conductive  up  to  24,  000°K  are  computed  This  upper  limit  was  decided 
upon  because  it  coincides  with  the  limit  of  Gilmore's  computations  of  the 
equilibrium  composition  of  air  (14). 

The  usual  practice  in  computing  the  conductivity  of  partially 
ionized  gases  is  to  add  the  Maxwell -averaged  total  electron  collision  cross 
sections  (weighted  with  the  number  density)  of  the  neutral  molecules  .o  a 
corresponding  equivalent  cross  section  of  the  ions,  the  latter  being  com¬ 
puted  on  the  basis  of  a  theory  which  applies  to  fully  ionized  gases.  The 
degree  of  approximation  involved  in  this  assumption  is  not  known.  In 
the  case  of  air.  .•»  further  element  of  uncertainty  is  introduced  by  the  fact 
that  the  cross  sections  of  some  of  the  neutral  species,  notably  N,  are 
not  known  with  certainty. 
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The  conductivity  cr  is  expressed  by 

2 

n  e 

o-  =  - — -  (A-l) 

m  C  T  n.Q. 
e  e  A-  j  j 

j 

where  is  the  number  density  of  free  electrons,  e  the  electronic  charge, 

m  the  electronic  mass,  n.  the  number  density  of  the  species  j  ,  Q.  its 

J  J 

Maxwell -averaged  total  electron  collision  cross  section,  and 


(k  =  Boltzmann's  constant,  T  =  temperature)  the  mean  speed  of  the  elec¬ 
trons.  Thermodynamic  equilibrium  with  an  electron  temperature  equal  to 
the  gas  temperature  is  assumed. 

The  summation  in  Eq.  (A-l)  is  extended  over  those  species, 
neutrals  as  well  as  ions,  which  contribute  appreciably  to  the  total  cross 
section.  An  inspection  of  the  equilibrium  composition  of  air  in  the  range 

from  6,  000°K  to  24,  000°K  and  a  density  ratio  p/  pQ  (pQ  =  standard  den- 
sity)  varying  from  10  to  10  indicates  that  the  neutrals  of  major  import¬ 
ance  are  N^.  N  and  0  (14).  In  view  of  the  approximate  nature  of  the 
computation,  the  inclusion  of  0^,  NO  and  other,  even  less  frequent  species, 
would  not  be  warranted.  The  cross  sections  for  N^»  N  and  0  above  an 
electron  tempera  .ure  of  6,000°K  are  believed  to  be  quite  similar  and  nearly 
independent  of  the  energy  in  the  range  considered.  (At  much  higher  tem¬ 
peratures,  the  effect  of  the  neutrals  becomes  negligible).  The  computations 

can  be  based  therefore  on  a  common  cross  section  of  the  neutrals,  Q  , 

n 

-15  2 

which  was  taken  as  0. 80  x  10  cm  .  In  view  of  the  uncertainty  in  the 
cross  sectional  data  of  N  and  0,  a  more  detailed  computation  would  not  be 
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justified. 

The  sum  over  all  cross  sections,  in  Eq.  (A-l),  can  therefore 
be  expressed  as 


T  n.Q.  =  n  Q  +  n.Q. 
v  j  j  n  n  li 


=  n 


A,  o 


£ 

Po 


(N,  +  N  +  0)Q  +  e'  Q. 

c  n  1 


(A-2) 


where  and  n.  are  the  number  densities  of  neutrals  and  ions,  respectively. 

n.  is  the  number  density  of  “air  atoms’*  at  standard  condition;  i.e. 

A  i  O 

19-3 

n.  =  5.38  x  10  cm  .  N,,  N,  0  and  e  are  the  number  of  N-mole- 

A  •  o  m 

cules,  N-atoms,  0-atoms,  and  free  electrons,  respectively,  per  “air 
atom”.  The  values  computed  by  Gilmore  have  been  used  for  these  fracti 

Within  the  limits  of  temperature  and  density  considered  1  re, 
the  effect  of  doubly  ionized  ions  on  the  conductivity  amounts  to  less  than  8% 
and  is  neglected.  Similarly,  the  effect  of  0  is  negligible.  Consequently, 
for  a  neutral  plasma. 


n. 

l 


n 

e 


n 


A,  o 


e 


(A-3) 


a  result,  which  was  already  utilized  in  the  derivation  of  Eq.  (A-2). 

Spitzer  and  Harm’s  (15)  results  for  the  conductivity  of  a 
fully  ionized  gas,  if  written  in  terms  of  an  equivalent  cross  section  Q.  for 
the  positive  ions,  give 


Si  ■ 


ire 

WT 


2  2  ln(<,Ce2) 


(A- 4) 
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for  singly  ionized  ions,  where 


m 


q  = 


2e 


kT 

2irn 


1/2 


e  J 


and  y_.  =  0.  582.  C  is  the  rms  electron  velocity, 
\b  e 


C  =  (— ) 

e  m 

e 


1/2 


interval,  the  computed 
Lamb  and  Lin  (13).  The 
ter.  However,  this  agree- 
es  were  used,  which  are 


The  results  of  a  computation  of  cr,  based  on  Eqs.  (A-l),  to 
(A-4)  are  given  in  Fig.  9.  It  is  noted  that  the  dependence  upon  density  is 
relatively  weak,  and  that  its  gradient  reverse s  withir.  the  considered  interval 
of  temperature.  Above  24,  000°K  the  effect  of  th <  econd  ionization  becomes 
increasingly  important. 

At  the  lower  end  of  the  tempera 
conductivity  agrees  with  the  experimental  dat 
computed  values  fall  inside  the  experimental 
ment  in  to  be  expected,  since  cross  sectional 
derived  from  these  data. 

Gilmore's  results  for  the  thermc  nic  variables  of  state 
of  the  air  behind  a  normal  shock  resulting  from  h>h  sonic  velocities,  have 
been  used  to  calculate  the  conductivity  for  different  flight  velocities  V  and 
altitudes  H  above  sea  level  (Fig.  10).  For  comparison,  the  velocity  of  a 
satellite  in  a  circular  orbit  of  7  x  10^  m  radius,  and  the  escape  velocity, 
have  also  been  indicated. 

In  many  applications  of  magnetoaerodynamics,  the  parameter 
crB^l/ pu  plays  an  important  role,  since  it  is  a  measure  for  the  ratio  of  the 
magnetic  forces  to  the  inertial  forces.  a,  p  and  u  are  defined  as  the  con¬ 
ductivity,  density  and  velocity  behind  a  normal  shock.  B  is  the  magnetic 
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induction,  and  1  a  characteristic  dimension  associated  with  the  body  pro¬ 
ducing  the  shock  wave.  In  general,  this  parameter  must  be  at  least  of 
order  one,  if  the  magnetic  field  is  to  produce  an  appreciable  effect. 

In  order  to  permit  a  quick  estimate  of  the  typical  field  strength 
required  in  magnetoaerodynamics  (assuming  thermodynamic  equilibrium). 
Fig.  11  has  been  drawn.  In  this  figure,  the  required  field  strength  is  given 
in  the  case  where  the  value  of  this  parameter  is  1,  and  for  a  length  1=1  m. 
Since  from  conservation  of  mass 

r 

p“  =  p^v 

where  p ^  =  p^(H)  is  the  free -stream  density,  the  computation  of  B  as 
a  function  of  H  and  V  is  straight-forward. 

Eq.  (A-l)  gives  or  correctly  only  in  the  case  of  weak  fields. 
Otherwise,  the  conductivity  is  no  longer  described  by  a  scalar  quantity. 

This  occurs  if  the  cyclotron  frequency  w  of  the  electrons  is  comparable  or 
larger  than  their  collision  frequency  v*  .  From 

Be 

u  =  — 
e  m 

e 


it  follows  from  equating  u  and  V  ,  that 

n  e 

B*  =  -S- 

9 


(A- 5) 


where  B'  can  be  considered  as  a  rough  estimate  of  the  maximum  field 
strength  B  for  which  Eq.  (A-l)  i*  still  valid.  The  locus  of  points  in  the 
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H  -  V  diagram  for  which  this  limiting  condition  exists,  is  indicated  in  the 
figure.  To  the  left  of  this  curve,  w  »-  V  . 
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Captions  to  Figures 


Fig. 

1 

The  functions  F(^),  G(  ^  )  and  their  derivatives,  for 

=  (3  =  1 

Fig. 

2  : 

Function  F^(^) 

Fig. 

3  : 

Function  F^'  {  £ 

Fig. 

4  : 

Function  F^(  £ 

Fig. 

5  : 

y 

Function  F^'  (  S 

Fig. 

6  : 

Function  G^{  £  j) 

Fig. 

7  : 

Magnetic  cushion  thickness  £ 

Fig. 

8  : 

Nusselt  numbers  Nu^  and  Nu^,  for  Pr  =  0.70 

Fig. 

9  : 

Estimated  electrical  conductivity  of  air 

Fig. 

10  : 

Estimated  electrical  conductivity  of  air  behind  a  normal 

shock 

Fig. 

11  : 

2 

Required  magnitude  of  magnetic  field  for  aB  1/pu  =1 

and 

1  =  1  m,  behind  a  normal  shock  in  air. 

Fig. 

12  : 

Model  used  in  glow-discharge  experiments  in  a  supersonic  stream 

Fig. 

13  : 

Typical  current- voltage  characteristic  for  glow- discharge  in 

supersonic  air  stream.  Copper  electrodes,  parallel  to  free 
stream  velocity.  Cathode  area  =  3cm^. 
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